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I 1. PROBLEMS; No problems were encountered during our analyses, 

2. ACCOMPLISHMENTS; 

1, Considerable effort was expended during i;he interval from 
11/1/81 to 4/30/02 in revising the manuscript that we submitted to 
the Journal of Geophysical Research in November, 1981. Although the 
paper was accepted 1n its original fprm, we felt that revision was 
needed. The manuscript deait with the discovery of a Precanibrian 
rift running NW-SE through Missouri w'is seen in both free air and 
Bouguer gravity anomalies and HCMM data. The revised manuscript is 
attached. Ue have substantially changed the paper, including the 
addition of magnetic field anomalies, basement rock types, and 
seismic refraction profile data for the area. Of importance to this 
study, we substantially revised the section on 1 inears seen in 
processed HCMM data covering the area. Included in the revised 
paper is a discussion of what HCMM data were utilized and a map that 
shows linear features seen from a merged daytime thermal image with 
a shaded relief map of topography. Also drawn on the Figure (Figure 
11 of the paper) are mapped faults and folds, and the outline of the 
rift running through the state. HCMM linenrs were mapped separately 
by three geologists, with the 1 inears drawn on the figure being ones 
seen by all three people, HCMM linears in some cases correspond to 
mapped features, in some cases they are extensions of mapped 
features, and in some cases they do not correspond to any mapped 
structure. The linears correlate nicely with the trend of the 
Precambrian rift as seen in the gravity data, suggesting that the 
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Unears ara Indeed related to the underlying basement structure. 

The pattern of faulting is consistent with uplift, suggesting that 
the rift is isostatically readjusting. HCMM data proved to be 
particularly useful because the synoptic coverage and the broad 
pixel width In effect filtered out all but regional-scale 
structures. 

The reason that JGR paper took so long to revise and the reason 
that these quarterly reports are delayed is that we were invited by 
the editor of EOS (Transactions, American Geophysical Union) to 
submit a paper on our Missouri rift, stressing digital image 
processing ofi potential field and topographic data. That paper is 
also attached - it appeared in the May 4, 1982 issue of E0S» with 
our color gravity maps on the cover. 

2. Attached are HCMM daytime thermal IR, night time thermal 
IR, and apparent thermal inertia images that have been 
contrast-enhanced and transformed to Mercator projections, Also 
attached is a daytime IR overlain onto a shaded relief map and an 
apparent thermal inertia image overlain onto a shaded relief map. 
Frame numbers are A-A0045-19420-2, A-A0045-19420-1, and 
A-A0044-08310-3 . These products were processed and used as part of 
a Masters Thesis finished in May 1, 1982 by dohn Strebeck. In the 
thesis Strebeck examined correlations between the HCMM data 
products, linears, and geologic units. The significant results are 
included in the attached J.6.R. paper. The thesis reference would 
be: 

Strebeck, O.W., 1982, Structure of the Precambrian basement in the 
Ozark Plateau as inferred from gravity and remote sensing data. 


Master of Arts Thesis, Washington Univ, , p., avaHable 

Univ. Microfilms, Ann Arbor, Michigan. 
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Strebeck now works with Mobil Oil Corp., Dallas, Texas, where he 
should be joining their remote sensing/ image processing group after 
preliminary “indoctrination" into the company. 

3. We have begun to examine the difference in information 
content between day IR, night IR, albedo, and thermal inertia 
images. For the frames listed above, the day IR image was found to 
have the greatest discriminability for rock types and with regard to 
structural patterns. The reason seems to be primarily due to 
topography in that most geologic units in southern Missouri exhibit 


a distinctive topography and most structural features that we have 
mapped are controlled by the distribution oF stream valleys. The 
day IR image best expresses these distributions because the 
brightness (relative temperature) is sensitive to differential solar 
heating due to slopes. For the three frames discussed above, the 
day IR contained 0.63 of the fractional variance of the three data 
sets, while the night IR contained 0.12 and the albedo data 
contained 0.25. Apparent thermal inertia images should provide an 
indication both of topographic roughness (since slope effects on 
differential solar heating have not been removed) and an indication 
of thermal anomalies that are controlled by variations in soil or 
rock exposure, vegetation type, and moisture content, We are 
presently concentrating on delineation of these latter features, 
i.e., identifying linears in southern Missouri that are not directly 
controlled by topography. 
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4. We have also coded the thermal IR and albedo data as hue, 
saturation, and brightness values to generate a color display, 
following the work of Haydn, University of Munich. However, we find 
that the transformation used by Haydn to go from hue, saturation, 
and brightness to red, green, and blue images is only correct for 
relatively narrow hue range. A universal derivation is attached, 
along with a color print where the day, night, and albedo images 
have been encoded as hue, brightness, and saturation. A color print 

is also included that has the apparent thermal inertia coded in a 
color spectrum and overlain onto a shaded relief map. 

3. FUTURE WORK: 

1. Presentation of results of analyses of HCMM data as part of a 
paper to be presented at COSPAR, Ottawa, Canada, on May 19, 

1982, entitled: 

Arvidson, R.E., E.A. Guinness, O.W. Strebeck, 1982, Structure of 
the Midcontinent basement - Topography, gravity, seismic, and 
remote sensing data. 

and included in the session on Geophysical Measurements of Major 
Crustal Features From Space. 

2. An oral presentation of our HCMM work at the spring meeting of 
the American Geophysical Union, Wash., DC. during the week of 
May 31, 1982. The abstract is attached and entitled: 

Strebeck, J.W., E.A. Guinness, R.E. Arvidson, 1982, Linear 
Structures in Missouri seen with Heat Capacity Mapping Mission 
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ABSTRACT 

Digital free air and Bouguer gravity anomaly images have been 
constructed for the region bordered by 25* to 49* N, lat, and 80’ to 
110* K. long, from approximately 287,000 station readings. The 
technique used to interpolate between station locations was based on 
a two dimensional spatial filter, where the average of the anomaly 
values located within the filter area was computed. The images 
contain as many as 256 contours (values in byte variable), so that 
subtle anomaly patterns can be identified and traced with much 
greater certainty than on most contour maps, A newly discovered 
feature in the midcontinent is a gravity low that begins at a break 
In the midcontinent gravity high in S.H. Nebraska, extends across 
Missouri in a NW-SE direction, and intersects the Mississippi Valley 
gtrbw to form the Pasco! a arch. The anomaly varies from 120 to 160 
km in width, extends about 700 km, and has a medial gravity high for 
, part of its length. The maximum Bouguer amplitude of the anomaly is 
approximately 34 milligals below values for the surrounding region 
and cannot be explained on the basis of a thickened section of 
Paleozoic sedimentary rock. The gravity data and the sparse seismic 
refraction data for the region are consistent with an increased 
crustal thickness beneath the gravity low. Discrete positive 
magnetic anomalie$, some of which correspond to mafic intrusives, 
are concentrated along the borders of the gravity low. Digitally 
enhanced thermal infrared images from the Heat Capacity Mapping 
Mission show a distinct alignment of linear structures with the 
gravity feature. The 1 inears in some cases correspond to mapped 
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high angle normal faults, to drape folds over relief within the 
Precambrian basement, and in some cases to extensions of mapped 
structures. The gravity anomaly also cuts across the major 
Precambrian boundary in S.E. Missouri marking the change from older, 
sheared granites and metasedimentary rocks, to younger granites and 
rhyolites. Given the cumulative evidence, the gravity anomaly is 
probably the present expression of a failed arm of a rifting event, 
perhaps one associated with the spread ’ng that led to or proceeded 
formation of the granite and rhyolite terrain of southern Missouri. 

INTRODUCTION 

We have been pursuing relationships between the pattern, age, 
and origin of structural features within the Precambrian basement 
rocks of southern Missouri and the locations and genesis of ore 
deposits. The area is a well kno'^ lead mining district, where 
Mississippi Valley type pb-Zn-Cu ores have accumulated in Cambrian 
carbonate rocks associated with stromatol itic reef and backreef 
facies (Gerdemann and Meyers, 1972). The location of these facies 
was controlled by the location of the shore line, which was in turn 
controlled by the pattern of faulting of the Precambrian basement 
(Grundmann, 1977; Sweeny et al., 1977; Evans, 1977; Paarlberg and 
Evans, 1977; Mouat and Clendenin, 1977). In addition, iron ores of 
magmatic origin can be found in the Precambrian basement rocks along 
fracture zones (Kisvarsanyi , 1976). 

Structural studies of basement rocks in southern Missouri have 
been pursued for a considerable amount of time (see: Kisvarsanyi 


PAGE 4 


and Kisvarsanyi, 1976), and a significant amount of information has 
been gained on the distribution of Precambrian rock types and ages 
(Kisvarsanyi , 1974; Bickford et al., 1981; VanSchmus and Bickford, 
1981). However, little has been done in terms of understanding how 
the structure of the region is related to the overall structural 
configuration of the midcontinent or how such features are related 
to the structural history of the area. In this paper we utilize 
digital image processing techniques to reduce and display a variety 
of potential field, topography, geologic, and remote sensing data 
for the midcontinent. The intent is to delineate how structural 
patterns in southern Missouri are related to the broader features of 
the midcontinent. 

DISCRIPTION OF DATA AND PROCESSING METHODS 


Digital image processing techniques have a potentially wide 
range of utility for display and analysis of geographically (i.e,, 
array) oriented data sets. In our case, data covf>ring the 
midcontinent were processed on a PDP-11/34 minicomputer with 
interactive image display peripherals, using standard digital image 
enhancement, filtering, and geometric operations. The reader is 
referred to standard texts such as Moik (1930) for further 
information on image processing techniques, and to Arvidson et 
al . (1982) for examples of the utility of image processing in 
processing, display, and interpretation of topographic and gravity 
data for the continental United States. 

National Oceanic and Atmospheric Administration (NOAA) digital 
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topography, together with land station measurements of gravitational 

acceleration, comprise two important data sets that we employed in 

the study. The area chosen for analysis of the topography and 

gravity ranges from 25° to 49“ N, lat. and 80“ to 110“ W. long,, 

thereby covering regions from the eastern edge of the Rockies on the 

west to the Appalachians on the east, and from the Superior Province 

on the north, to the Gulf of Mexico on the south. The topography 

» 

consists of average elevations for areas covering 30 seconds in both 
latitude and in longitude. Gravity stations are typically spaced 
three km apart, but vary between hundreds of meters to ten 
kilometers over the study area. The gravity data were reduced to 
free air and to Bouguer anomalies courtesy of the Defense Mapping 
Aerospace Agency Center, St. Louis, Missouri. The reference field 
at sea level used was based on the 1967 International Gravity 
Fonnula, with all the data referenced to the 1971 International 
Gravity Standardization Net. Bouguer anomalies were computed based 
on a slab model with a density of 2.67 gm/crn3. No local terrain 
corrections were included. 

Magnetic anomalies were. also included in the analysis. 
Unfortunately, digital data for the region were not available at a 
reasonable cost. Consequently we restricted inclusion of magnetic 
coverage to the state of Missouri, based on the 1943 statewide 
contour map of vertical field intensity anomalies. The state map 
was photographed onto film and then the film was digitized, 
resulting in a digital image of magnetic anomaly contours for 
Missouri. Similarly, the basement rock type map for the 
midcontinent compiipd by Bickford et al. (1981) and VanSchmus and' 
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Bickford (1981) was transformed to a digital format for use in 
comparison with other data sets. 

Other digital data sets that were utilized in this work 
included Heat Capacity Mapping Mission (HCMM) images. HCMM 
consisted of a satellite with sensors capable of imaging the surface 
in the visible to the reflected infrared (0.5 to 1.1 micrometers) 
and in the thermal infrared (10.5 to 12.5 micrometers) (Price, 

1977). Each image element in the HCMM data covers about 500m across 
and one scene covers about 700 km in width and in height. The 
visible-reflected infrared sensor was used to measure the broad band 
albedo of the surface, while thermal infrared data acquired during 
the day and night provided information on the magnitude of diurnal 
temperature changes. These parameters can be used to solve for the 
thermal inertia of the surface, a thermophysical property that 
depends on the density, thermal conductivity, and specific heat 
(Kahle et al., 1976). Thermal infrared data, including estimates of 
thermal inertia, have been shown to be useful in delineating 
structural features in a variety of geological contexts (Sabins, 
1969; Offield et al., 1975; Watson, 1981). HCMM data were 
included in this study both because of the applicability of thermal 
data in structural studies and because the synoptic view is 
appropriate for examining the surface expression of regional-scale 
structures. HCMM images for southern Missouri were contrast 
enhanced, digitally registered to one another, and displayed in a 
variety of formats designed to emphasize structural trends. Details 
of HCMM processing techniques will be discussed in a later section. 

The first step in the analysis of topographic and gravity data 
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was to scale the dynamic range of the data to fit within the range 
of a byte variable (i.e., 8 bits or 256 discrete values). The byte 
data were then stored in image arrays, registered to a common base 
if needed, and processed. Final displays were transformed to 
Mercator projections and overlayed with latitude-longitude grids. 

The byte-encoded topographic data were first stored in an array with 
30 second spacing in both latitude and longitude. The data could be 
displayed as an image if the byte-encoded values were converted to 
brightness or to color values. Alternatively, a shaded relief image 
could be constructed if a photometric function for the surface is 
assumed. Figure 1 shows a shaded relief image transformed to a 
Mercator projection and depicting topography under the assumption 
that the surface obeys Lommel-Seel iger scattering law (Batson et 
al., 1975). The simulated sum is from the northeast at 20" above 
the horizon. 

Generation of images from the gravity anomaly values is more 
complicated than generating displays from the topographic data. The 
reason is that the gravity stations are not located along regular 
grid intersections. We chose a simple, but pow,'-ful spatial 
filtering technique for interpolating between station locations. 

The byte-encoded anomaly values were first assigned to array 
locations closest to the station locations. This procedure produced 
an array in part occupied by valid data and in part occupied by 
blank zones. The filtering algorithm that we then applied was 
developed by Eliason and Soderblom (1977), used to produce 
topographic maps from the Pioneer-Venus altimetry data (Pettengill 
et al., 1980), and applied to generating gravity images for the 
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continental United States by Arvldson et al. (1982), The technique 
involves use of a spatial filter of N x N elements. The average of 
the valid data points within any given filter location is computed 
and used to replace the midpoint value, but only if a valid datum 
does not already exist at that array location. In practice, the 
gravity data were processed using several filter passes, beginning 
with a 3 X 3 element filter and ending with a 21 x 21 element 
filter. The choice of filter sizes’ was governed by station 
spacings, which varied from hundreds of meters to ten kilometers. 

The result Is an interpolated data set, except for regions that had 
so few valid data points that a user-defined threshold was not met. 
In our case, we set the threshold so that at least 20 % of the 
elements for any filter position would have to have been occupied by 
valid data for an interpolation calculation to proceed. 

The interpolated gravity data can be displayed as a gray tone 
image, as is done in Figure 2 for free air anomalies. Black areas 
in the image correspond to zones with too few stations to allow 
valid interpolations. The gravity data can also be displayed as a 
shaded relief image, where gravity highs and lows are illuminated as 
If they were hills and valleys. Figure 3 is a shaded relief image 
of the free air anomalies with the simulated sun placed at 15“ above 
the northeastern horizon. 

The gravity anomaly images displayed in Figures 2 and 3 and 
equivalent Bouguer anomaly presentations have been checked against 
published maps (e.g. Woolard and Joesting, 1964; McGinnis et al., 
1979; Simpson and Goodson, 1981). Generally the images and 
published maps correspond. However, the images have intrinsically 
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more information displayed, for at least three reasons. First, the 
grid spacing used to interpolate between station readings is more 
closely spaced than the spacing used to generate most published 
continental-wide contour maps. Second, there can be as many contour 
intervals in our displays are there are values in a byte (256). 
Third, because of the particular filtering algorithm used, areas 
with closely spaced stations retain local details of anomaly 
patterns. For instance. Figures 2 and 3 show a major gravity 
feature (free air low) that begins at a break in the midcontinent 
gravity high and extends about 700 km to the southeast. One of the 
distinguishing aspects of this feature is not the amplitude of the 
low but rather the sharp gravity gradient associated with the edges 
of the anomaly. Although sections of the feature have been noted in 
the past (see: Phelan, 1969; Cordell, 1979; Russ, 1981), the 
images shown in Figures 2 and 3 provide new information and a 
broader perspective that show the trend and location much more 
accurately than can be seen in published maps. 

Finally, an important aspect of comparing data sets is the 
ability to merge or overlay one data set onto another to produce a 
visual display that maximizes the eye's ability to see correlations 
between the data sets. A useful technique is to allow one data set 
(gravity anomalies, for instance) to control the hue (dominant 
wavelength) and saturation (degree of purity) of a color image, 
while another data set (topography, for instance) controls the color 
brightness (Arvidson et al., 1982). Also, it is possible to combine 
values of a given parameter with a shaded relief presentation of 
that parameter by allowing the value to control the color hue and 
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saturation, while the local gradient, expressed in shaded relief 
form, controls the brightness (Pettengill et al., 1980; Kobrik, 
1982). For example, Figure 4 is a version of the free air image 
where the anomaly values have been color-coded and overlain onto a 
shaded relief version of the anomalies. The effect is to produce an 
enhancement showing information related both to the value of the 
anomaly and to the local gravity gradient. 

BASEMENT STRUCTURE IN MISSOURI AS DEFINED 

BY GRAVITY ANOMALIES AND SEISMIC PROFILES 

Figure 5 is a sketch map showing the major structural features 
seen in the free air anomaly images shown in Figures 2, 3, and 4 and 
in Bouguer images covering the same region. The free air data 
facilitated mapping of major topographic features associated with 
the Cordillera, Ouachitas, and the Appalachians. The Bouguer data 
provided new information that allowed us, for instance, to map the 
strong positive anomalies associated with the Mississippi embayment. 

The midcontinent gravity high is a major gravity anomaly in the 
midcontinent. In our data this feature has a maximum free air 
anomaly of 85 milligals and a maximum Bouguer anomaly of 50 
milligals. The midcontinent gravity high is about 70 km wide, 
contains Keewanawan basalts that date at 1.1 Ga, and has been 
interpreted as a failed continental rift (Chase and Gilmer, 1973). 
Flanking lows on either side of the high have been modeled as thick 
(several kilometers) arkosic sediments deposited as the load 
associated with the basalts caused regional subsidence. 
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Considerable structure an also be seen in the greenstone-granitic 
terrain of the Superior Province located to the northwest of the 
midcontinent gravity high. The Wisconsin Arch, the Ouachitas, and 
the Wichita-Arbuckle system are just a few of the other features 
that have a recognizable gravity signatures. For reference, the 
Quachitas have a maximum free air and Bouguer anomaly magnitudes of 
76 and 54 milligals below the surrounding regions. 

A subtle, but pervasive gravity feature seen in both the free 
air and Bouguer images is a linear low that begins at a break near 
the southeastern end of the midcontinent gravity high (40,5“ 

N. lat., 96“ W. long.), strikes in a southeasterly direction, and 
extends to the Mississippi Embayment (36“ N. lat., 90“ W. long.). 

The low varies between 120 to about 160 km in width, extends about 
700 km in length, and for part of its length exhibits a medial 
gravity high. The maximum Bouguer anomaly associated with the low 
is about 34 milligals below values for surrounding regions. The 
feature begins just to the northwest of Missouri and ends just 
beyond the southeastern boundary of the state. We therefore 
informally refer to the feature ar the Missouri gravity low. As 
discussed in Arvidson et al. (1982), the intersection of the 
Missouri gravity low with the Mississippi Valley graben as defined 
by Kane et al . (1981), is the site of the majority of the 
microseismic epicenters recorded in the 1970’ s by the St. Louis 
University seismic network (Stauder et al . , 1977). In addition, the 
northern boundary of the intersection is the site of the circular 
Bloomfield anomaly, while the southern boundary is the site of the 
Covington anomaly. Both features can be seen in the free air and 
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Bouguer data as distinct, circular positive anomalies. Kane et 
al. (1981) suggest that these two anomalies, which are also positive 
magnetic anomalies, are due to intrusions of magmas during the 
Mesozoic Era. The intersection of the Missouri gravity low and the 
Mississippi Valley graben is also the site of the Pascola arch as 
defined by Phelan (1969) and Ervin and McGinnis (1975). 

Several gravity profiles were generated across the Missouri 
gravity low to illustrate its form and to model subsurface density 
configurations. The location of the southern-most profile is shown 
in the sketch map in Figure 5. Values for the topography, free air, 
and Bouguer anomalies for this profile are shown in Figure 6. A 
simple model for the anomaly would be a thickened section of 
relatively low density sedimentary rock overlying a downwarped 
region of the Precambrian basement. However, drill holes covering 
the general area of the profile A-A' indicate a typical cover of 
Paleozoic sediments of only a few hundred meters over the low, with 
no discernable thinning on either side (Kisvarsanyi , 1974). Thus, 
the gravity signature must be related to a inhomogeneity within the 
Precambrian basement r>cks, Cordell (1979) reached a similar 
conclusion for that part of the low located in the southwestern part 
of the Rolla quadrangle (37" to 38" N. lat.; 90" to 92" W. long,), 
v^here there is a slight local thickening of sedimentary cover over 
the gravity low. 

Some seismic data exist for Missouri that provide information 
on the possible subsurface configurations of the crust that would 
give rise to the observed gravity anomaly patterns. Stewart (1968) 
conducted a reversed seismic refraction survey of about 300 km in 
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length along an east-west line in northern Missouri crossing the 
gravity low at about 39.8* N. lat. Results Indicate that the crust 
consists of three major layers, with average depths of about 5, 20, 
and 40 kin and a 3/4 degree of dip toward the west, Stewart (1968) 
also conducted a reversed profile in southern Missouri that extends 
across the gravity low at about 50 kin to the northwest of the 
gravity profile A-A* (Figure 5), The refraction data for Stewart's 
(1968) southern profile exhibit very low signal/noise ratios, 
perhaps because of significant lateral inhomogeneities within the 
crust. Results are compromised by the low quality of the data, 
although an arguement can be made for a thickened crustal section in 
the southwestern half of the profile and somewhat higher crustal 
velocities beneath the northeastern half of the traverse (Nuttli, 
I 9 V 6 ). Finally, McCamy and Meyer (1966) conducted a seismic 
refraction survey from Little Rock, Arkansas to Cape Girardeau, 
Missouri, reversing a survey done several years earlier. The 
refraction profile cut across the Bloomfield gravity anomaly and ran 
in a NE-SW direction within the Mississippi ernbayinent (i.e., to the 
southeast of profile A-A* of Figure 5). As discussed by Ervin and 
McGinnis (1975), this region is unusual in that the crust is 
slightly thicker and underlain by an anomalously high velocity layer 
as compared to Stewart's (1968) profile in northern Missouri. The 
anomalous layer probably corresponds to high density material at the 
base of the crust, material that may have been emplaced in 
association with the Precainbrian rifting that produced the Reelfoot 
rift. 

In summary, the seismic data suggest a slightly thicker crustal 
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section under the Missouri gravity low, although this interpretation 
is compromised by poor data and the complicating Influence of rocks 
associated with processes that probably lead to formation of the 
Mississippi embayment. If it is assumed that the Missouri gravity 
low is due to a thickened crust, then a 34 milligal anomaly, and a 
density contrast of 0.3 gm/cm3 between crust and mantle, is 
consistent with a crustal excess of about 3.3 km under the anomaly 
(Strebeck, 1983). Alternatively, if it is assumed that lateral 
density variations cause the gravity anomaly, than the anomaly would 
be consistent with rocks that were 0.1 gm/cm3 less dense than 
surrounding materials for the first 4 to 8 km below the surface 
(Strebeck, 1982). Of course, other models of crustal 
inhomogeneities can be generated. However, in the absence of other 
constraints, they are not worth pursuing at this point. 

j 

I 

RELATIONSHIP OF MISSOURI GRAVITY LOW TO I 

PRECAMBRIAN ROCK PROVINCES AND TO ! 

i 

MAGNETIC ANOMALIES 

Information on the distribution of Precambrian rock types in 
Missouri is limited because of the relatively few drill holes that 
have penetrated into the basement and because of the complexity of 

, 

the Precambrian geology in the area (Kisvarsanyi , 1974,* Bickford et 

al., 1981; VanSchmus and Bickford, 1981). There is, however, a 

major boundary between 1.6 Ga old rocks composed largely of sheared 

granites and metasediments, and younger 1.4 6a old granites and 

rhyolites. The boundary runs in a NE-SW direction, as is shown in j 
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Figure 7, where the basement rock map of Bickford et al. (1981) has 
been overlain onto a shaded relief map depicting Bouguer anonalles. 
The Missouri gravity low cuts across the boundary at nearly right 
angles. There 1s also clearly an extension of what are thought to 
be older metasedimentary rocks Into the younger terrain along the 
gravity low. Kisvarsanyl (1974) notes that the area containing the 
reentrant of older rocks is also a structural high on contour maps 
depicting the paleotopography of the Precambrian surface and on 
structured contour maps of Paleozoic sedimentary formations. 
Presumably, older Precambrian rocks are exposed along the low 
because of Isostatic readjustment and erosion of the younger 
granites and rhyolites preferentially along the feature. 

The observation that the gravity low cuts across the major 
Precambrian age and r-ocK type boundary in Missouri suggests that the 
low is related to a structural feature. Further support for this 
suggestion can be found in the locations of areas with discrete, 
positive magnetic anomalies. Figure 8 shows regions with vertical 
intensity magnetic anomalies higher than 600 gammas as black 
splotches superimposed on the gray tone and shaded relief versions 
of Bouguer anomalies for Missouri. The 600 gamma contour interval 
was chosen on the basis of delineating discrete magnetic highs. The 
magnetic highs are in part clustered along the flanks of the gravity 
low, i.e., where the gravity gradient is high. Other magnetic highs 
coincide with discrete gravity highs, such as in the area underlain 
by the Bloomfield intrusion in southeastern Missouri or the Salem 
Intrusion in the middle of the gravity low (Figure 8), Clearly, 
more work needs to be done to establish the pattern of magnetic 
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anomalies in more detail, along with establishing the relationship 
of the anomaly patterns and the gravity low. Examination of 
unpublished maps constructed by I, Zietz supports the observation 
that the magnetic highs are along the flanxs of the gravity low. 

Previous work shows that some magnetic highs in Missouri can be 
correlated with relatively iron rich igneous rock bodies (Phelan, 
1969; Cordell, 1979). The location of the discrete magnetic highs 
long the flanks of the gravity low suggests that the flanks 
correspond to zones of weakness where magmas have been intruded, 

Kane et al. (1981) note a similar correlation of magnetic highs 
along tre perimeter of the Mississippi Valley graben. The gravity 
and magnetic anomaly patterns, and the observation that the gravity 
low cuts across rock and age provinces, provide strong evidence that 
the Missouri gravity low is a consequence of a major inhomogeneity 
within the crust. 

RELATIONSHIPS BETWEEN FOLDS, FAULTS, AND LINEARS 
IN SOUTHERN MISSOURI AND BAi^cMENT STRUCTURE 

In this section we describe relationships between the Missouri 
gravity low and the patterns of faults, folds, and 1 inears that 
exist within the Paleozoic sedimentary rocks covering Missouri, The 
particular tack that we take is to map the patterns of linear 
features from HCMM thermal images, compare those 1 inears to the 
distribution of known faults and folds, and finally we relate the 
trends for the structural data to the gravity low. 

Figure 9 shows a HCMM thermal infrared image over southern 
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Missouri taken during relatively clear atmospheric conditions at 
approximately 1;30 p.m. on dune 10, 1978, The frame has been 
contrast-enhanced and transformed to a Mercator projection. Water 
appears dark (cool) because it does not warm as quickly as 
vegetation, soil, or rock. A correlation plot between the apparent 
temperature seen in this image and a map of the magnitude of 
topographic slopes demonstrates that the dominate control on 
temperature is relief. Areas with high relief show as relatively 
dark (cool), while relatively flat regions are bright (warm). For 
instance, the bright polygonal zone centered at 37.6“ N. lat., 90.5“ 
W. long, corresponds to flat terrain underlain in part by igneous 
rocks of the St. Francois Mountains and in part by the surrounding 
Cambrian sedimentary rocks. The dark area to the southwest of this 
region corresponds to Ordovician carbonates that have been dissected 
by streams. Flood plains of the streams and the major rivers are 
flat and therefore relatively bright (hot). Finally, the relatively 
smooth, bright regions to the north of the Missouri River and east 
of the Mississippi River coincide with the glacial till deposits 
related to the Wisconsin and earlier periods of glaciation. 

We also combined a shaded relief image of topography with the 
daytime thermal image by multiplying the two data sets on a element 
by element basis. The resulting product is shown in Figure 9. 
Clearly, linear features are enhanced in the merged data relative to 
either the shaded relief map or thermal images shown alone. The 
reasons are at least three fold. First, topographic linears can be 
readily discerned on shaded relief maps (Wise, 1976). Second, 
topographic linears associated with river valleys would be enhanced 
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because of the higher teiiiperatures for the valleys as opposed to the 
hills. Combining the thermal image with the shaded relief map 
provides, in effect, a nonlinear enhancement of the topography. 
Third, any linear thermal f(?atures not related to topography would 
remain in the data. This product was used as our primary base for 
constructing a 1 inears map for southern Missouri. An apparent 
thermal inertia image was also constructed using the formulation of 
Price (1977) and overlayed onto the .shaded relief image. However, 
this product did not provide a significant amount of additional 
information, proba»)ly because the night ime thermal and the albedo 
data contain only 25% of 13% of the variance inherent in the 3 
dimensional thermal day-thermal night-albedo data set. 

Figure 11 is a sketch map showing 1 inears from the processed 
HCMM data. Also shown are folds from a version of the Missouri 
structures map (McCracken 1971) and faults from the latest geologic 
map of the State (Anderson et al . , 1979). The folds are thought to 
be drape folds over the Precambrian surface in areas with 
significant relief due to basement faulting (McCracken et al., 

1971). Most of the faults are high angle normal faults, although 
some evidence for shear movement can be found (McCracken et al., 
1971). The area occupied by the high temperature region of the 
St. Francois Mountains and surroundings is also drawn on Figure 11. 
Note that linears, faults, and folds generally strike in a direction 
parallel to the gravity low. Also, the sharp southwestern edge of 
the St. Francois Mountains is coincident with the northern flank of 
the low. Further away from the low the the azimuths of features 
tend to disperse. Clearly, the basement in homogeneity associated 
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with the Missouri gravity low has had a pronounced influence on the 

pattern of faulting and folding that has propagated through the 

Paleozoic sedimentary cover. The preponderance oF norma’* faulting 

in the area implies, as does the reentrant of older basement rocks 

along the gravity low, that the dominate structural activity 

associated with the crustal inhomoqeneity beneath the gravity low 

has been one of vertical readjustment, preferential uplift along 

« 

the low is also consistent with the pattern of vertical 
displacements shown by the mapped faults in the area (Figure 11). 

THE AGE, ORIGIN, AND EVOLUTION OF THE 
MISSOURI GRAVITY LOW 

The Missouri gravity low must be older than the Paleozoic 
sediments that have been deposited over the reentrant of older rocks 
into the granite-rhyolite terrain. The reason is that any uplift 
and stripping of the younger granite-rhyolite rocks along the 
gravity low must have occurred before the Precambrian surface was 
buried. The oldest abundant sedimentary rocks in southern Missouri 
consist of Upper Canbrian sandstones called the Lamotte Formation 
(Anderson et al., 1979). 

A number of events affected the region now occupied by the 
gravity low during the Precambrian including; (1) The set of 
processes that led to formation of the granite-rhyolite terrain of 
southern Missouri. Bickford et al . (1981) speculate that a rifting 
event or a convergent plate margin might explain formation of the 
granite-rhyolite units; {?.) Rifting that was associated with the 
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opening of the Reel foot basin, the Precambrian precursor to the 
Mississippi enibayment (Ervin and McGinnis, 1975); (3) Rifting 
associated with formation of the niidcontinent gravity high (Chase 
and Gilmer, 1973); and (4) The postulated continent-continent 
collision associated with the Grenville orogeny (Burke and Dewey, 
1973). 

At this point it is useful to summarize the characteristics of 
the Missouri gravity lew. The feature has a negative Bouguer of 
about 34 mil'iigals, a width varying from 120 to 160 km, it has a 
medial high for part of its length, and it extends about 700 km in 
length. The low is perhaps best expressed within the younger 
granite-rhyol ite terrain of southern Missouri. The flanks of the 
feature exhibit magnetic highs, some of which are probably due to 
Igneous intrusions. Finally, the feature controls the pattern of 
normal faulting in the Precambrian basement rocks and overlying 
sedimentary cover. Except for the lack of a discernable thickness 
of sedimentary fill overlying the granite-rhyolite terrain, such a 
description is consistent with an origin as the failed arm of a 
triple junction. In fact, there are examples of relatively recent 
rifts that have negative anomalies with no associated sediment fill, 
although Burke and Whiteman (1973) interpret these anomalies as 
being due to ponding of magmas at the base of the crust. 

The Missouri gravity low could have formed as the failed arm of 
a triple junction that existed before formation of the 
granite-rhyolite terrain. The triple junction site would have been 
to the southeast, with the other two rifts becoming active spreading 
centers. Closure of that spreading center, with subsequent 
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continent-conbinerit collision, may have been the most plausible 
method of generating the abundant alkalic igneous activity that 
produced the younger granite-rhyolite terrain. In the process, most 
of the sediment fill within the failed rift would have been 
metamorphosed, intruded into, and covered by the alkalic magmas 
associated with formation of the granite-rhyolite terrain. The 
metasediments found in the reentrant of older rocks in this terrain 
(Figure 7) may be the remnants of this valley fill. Subsequent 
events, such as the formation of the Reelfoot rift that is 
postulated to have occurred at approximately 1.2 Ga (Ervin and 
McGinnis, 1975), could have served to reactivate parts of the rift. 
Reactivation certainly occurred during the mid-Paleozoic to 
Mesozoic, when the southeastern end of the low was structurally 
active as the Pascola arch. We stress that our interpretation , 
although plausible, is by no means unique and serves only as a 
working hypothesis to be updated or discarded as new seismic, 
potential field, and rock-type data are acquired. It is of interest 
to note that an origin as a failed rift arm may provide a new 
perspective on the source of the Pb-Zn-Cu in the Paleozoic 
sedimentary cover, since base metal mineralization is a common 
process associated with rifting events (Burke and Whiteman, 1973). 

SUMMARY AND IMPLICATIONS 

1. Standard digital image processing techniques have been used to 
interpolate and display free air and Bouguer anomalies and 
topography as gray tone, color-coded, and shaded relief images 
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for the miclcontlnent as part of an analysts of how structural 
features within the Precambrian basement of Missouri relate to 
the broader structure of the midcontinent. In addition, the 
Precambrian basement rock type map of Bickford et al . (1981) and 
magnetic anomalies for Missouri have been digitally merged and 
displayed with the gravity data to facilitate examination of 
correlations. Finally, Heat Capacity Mapping Mission (HCMM) 
thermal image data covering southern Missouri have been enhanced 
and overlayed onto shaded relief versions of topography. 

2. A major and previously unrecognized gravity anomaly (Missouri 
gravity low) in the midcontinent is a 700 km long, 120-160 km 
wide free air and Bouguer low, with a Bouguer amplitude of up to 
34 milligals. The anomaly displays a medial gravity high for 
part of its length. The low begins at a break in the 
midcontinent gravity high in southeastern Nebraska, extends in a 
southeast direction through Missouri, and intersects the 
Mississippi Valley graben as defined by Kane et al. (1981). The 
northern and southern edges of the intersection of the Missouri . 
gravity low and the Mississippi Valley graben are the sites of 
discrete positive gravity anomalies that are thought to be due 
to plutons. Also, many of the epicenters associated with the 
New Madrid seismic area are within the crustal block defined by 
the intersection of the gravity low and the Mississippi Valley 
graben. Discrete positive magnetic anomalies exist along the 
flanks of the Missouri gravity low. Finally, the gravity low is 
the site of a reentrant of older Precambrian metasedimentary 
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rocks into the younger granites and rhyolites that underlie much 
of southern Missouri, The amplitude of the low is too high to 
be related to thickening of Paleozoic sedimentary rocks. Rather 
the anomaly is probably due to a basement inhomogeneity. Two 
seismic refraction profiles ( Stewart » 1968) cut across the low 
in two places in Missouri, In both cases, the seismic data are 
consistent with a slightly thicker crustal section beneath the 
gravity low. The magnitude of the gravity anomaly is consistent 
with a crustal excess of 3.3 km as compared to surrounding areas 
or with a crust that is slightly less dense (0.1 gm/cm3 
contrast) in the uoper 4 to B km, 

3. A HCMM daytime thermal image was overlain onto a shaded relief 
image of topography covering the area underlain by the gravity 
low in southern Missouri. These data show that linears trend in 
the same direction as the low. On the other hand, linears tend 
to be dispersed in azimuth away from the gravity anomaly. 

Mapped faults (high angle normal faults, mainly) and drape folds 
over basement relief in the Paleozoic sedimentary cover of 
southern Missouri follow a similar pattern. Some of the HCMM 
linears correspond to the mapped structures, in some cases they 
are extensions of such features, and some linears do not 
correspond to mapped features. In any case, the correlation of 
the gravity low and the structural data suggest that the low is 
a major basement rift and that the structures in the Paleozoic 
deposits are due to vertical readjustments related to the 
crustal inhomogeneity that gives rise to the low. 
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4. The Missouri gravity low is older than the oldest Paleozoic 
sedimentary rocks that cover the reentrant of metasediineivtary 
rocks into the younger granites and rhyolites. The oldest rocks 
correspond to the Upper Cambrian Umotte Formation. The gravity 
low could be the present expression of a rift that formed as a 
failed arm of a triple junction occurring before or during the 
events that led to the younger granite-rhyolite terrain of 
southern Missouri. Emplacement of the igneous rocks would have 
effectively covered or metamorphosed sedimentary rocks that 
formed as valley fill with the rift. Reactivation of the rift 
has certainly taken place since it formed. Association of the 
Missouri gravity low with a unique history is difficult at this 
point because of lack of information on basement rock types, 
ages, and seismic controls on crustal configurations beneath the 
low. 


ACKNOWLEDGMENT 


This work was carried out with funding from the NASA Office 
of Space and Terrestrial Applications, Non-Renewable Earth 
Resources, through a contract (955959) from the Oet Propulsion 
Laboratory and by the Heat Capacity Mapping Mission Data 
Analysis Program, Goddard Space Flight Center under Contract 
NAS5-26533. Processing was conducted with a PDP'-ll/34 system 
located within the Washington University Regional Planetary 
Image Facility. Thanks are extended to St. Louis University for 
seismic data and to the Geopositional Division, Defense Mapping 
Aerospace Agency Center, St. Louis, Mo., for gravity data. 


PAGE 26 


REFERENCES CITED 

Anderson, K.H., 1979, Geologic map of Missouri, Missouri Geological 
Survey publication, 1 sheet. 

Arvidson, R.E., E.A. Guinness, 0. Strebeck, G. Davies, K.J. Schulz, 
1982, Image processing applied to gravity and topography data 
covering the continental United States, EOS, 63, 261-265. 

I 

Batson, R.M., K. Edwards, E.M. Eliason, 1975, Computer-generated 
shaded relief images, J. Res. U.S. Geol. Surv., 3, 401-408, 

Bickford, M.E., K.L. Narrower, W.J. Hoppe, B.K. Nelson, R.L. 
Nusbaum, tl,J, Thomas, 1931, Rb-Sr and U-Pb geochronology and 
distribution of rock types in the Precambrian basement of 
Missouri and Kansas, Geo. Soc. American Bull., 92, 323-341. 

Braile, L.W., W.J. Hinze, J.L. Sexton, 1981, Ancient rift complex 
and related seismicity in the northern Mississippi embayment 
area, Geol. Soc. Am, Abs. with Programs, 13, 415. 

Bretz, J.H., 1965, Geomorphic history of the Ozarks of Missouri, 
Missouri Geol. Survey and Water Resources, 2nd Series, 41, 134 p. 

Burke, K. and A.J. Whiteman, 1973, Uplift, rifting, and break-up of 
Africa, in D.H. Tarling and S.K. Runcorn, eds.. Implications of 
continental drift to the Earth Sciences, Academic Press, 735-755. 

Chase, C.6. and T.H. Gilmer, 1973, Precambrian plate tectonics: The 
midcontinent gravity high. Earth Planet. Sci. Lett., 21, 70-78. 




PAGE Z7 


Cordell, L,, 1979, Gravity and aeromagnetic anomalies over basement 
structure in the Rolla quadrangle and the S.E. Missouri lead 
district, Economic Geology, 67, 409-425. 

Dewey, J.F. and K. Burke, 1973, Tibetan, Varlscan, and Precambrlan 
basement reactivation; Products of continental collision, 0, 
Geology, 81, 683-692. 

Eliason, E. and L.A. Soderblom, 1977, An array processing system 
for lunar geochemical and geophysical data. Proceed. Lunar and 
Planet. Scl. Conf. 8th, 1163-1170. 

Ervin, C.P. and L.D. McGinnis, 1975, Reelfoot rift; Reactivlted 
precursor to the Mississippi embayment, Geol. Soc. An. Bull., 86, 
1287-1295. 

Evans, L.L., 1977, Geology of the Brushy Creek Mina, Viburnum 
Trend, Southeast Missouri, Economic Geology, 72, 381-390, 

Gerdemann, P.E. and H.E. Meyers, 1972, Relationships of carbonate 
facies patterns to ore distribution and to ore genesis In the 
S.E. Missouri lead district, Economic Geology, 67, 426-433. 

Grundmann, W.H., 1977, Geology of the Viburnum No. 27 Mine, 

Viburnum Trend, Southeast Missouri, Economic Geology, 72, 

349-364. 

Kahle, A.B. and A.R. Gillespie, and A.F.H. Goetz, 1976, Thermal 
inertia Imaging - A new geologic mapping tool, Geophys. Res. 
Lett., 3, 26-28. 


PAGE 28 


Kane, M.F., T.6. Hilclenbrand, and 0,0, Hendricks, X981, Model for 
the tectonic, evolution of the Mississippi einbaynent and its 
contemporay seismicity, Geology, 9, 563-568, 

Kisvarsanyl, E.6,, 1974, Operation basement; Buried Precambrlan 
rocks of Missouri - Their petrography and structure. Am, Assoc, 
Petroleum Geologists Bull., 58, 674-684, 

Kisvarsanyl, 6. and E.B. Kisvarsanyl, 1975, Ortho-polygonal 
tectonic patterns In the exposed and burled precambrlan basement 
of southeast Missouri, Proceedings of the First International 
Conf. on the New Basement Tectonics, Utah Geological Association, 
Pub. 5, 169-182. 

Kisvarsanyl, G., 1976, Precambrlan inetallogenesis in the St. 
Francois Mtns. igneous province, S.E. Missouri, in Kisvarsanyl, 
E.G., ed., Studies In Precambrlan geology of Missouri, Missouri 
Geological Survey, Rept. Inves. 61, 164-173. 

Kobrick, M., 1982, Topography of the terrestrial planets. 

Astronomy, 10, 18-22. 

McCamy K. and R.P. Meyer, 1966, Crustal results of fixed multiple 
shots in the Mississippi embayinent. Am. Geophys. Union Monograph 
10, 370-381. 

McCracken, M.H., 1971, Structural Features of Missouri, Missouri 
Geological Survey Rept. Inv., No. 44, 68p. 

McGinnis, L.D., M.G. Wolf, J.J. Kohsmann, C.P. Ervin, 1979, 

Regional free air anomalies and tectonic observations in the 


PAGE 29 


United States, J, Geophys. Res., 34, 001-601. 

Molk, J.6., 1980, Digital processing of remotely sensed images, 

NASA SP-431, 330 p. 

Mouat, M.M. and C.W. Clendenisri, 1977, Geology of the 0?:ark Lead 
Company Mine, Viburnum Trend, Southeast Missouri, Economic 
Geology, 72, 398-407. 

Nuttli, O.W., 1976, Seismological studies of Missouri crustal 
structure, in studies in Precambrian Geology of Missouri Div, 

Nat. Resources, Rept. Inves. 61, 184-189. 

Offield, T.W,, 1975, Thermal infrared images as a basis for 
structure mapping. Front range and adjacent plains in Colorado, 
Geol. Soc. Atn. Bull, 86, 495-502. 

Paarlberg, N.L. and L.L. Evans, 1977, Geology of the Fletcher Mine, 
Viburnum Trend, Southeast Missouri, Economic Geology, 72, 

391-397. 

Pettengill, G., E. Eliason, P. Ford, G. Loriot, H. Masursky, 1980, 
Pioneer-Venus radar results: Altimetry and surface processes, J. 
Geophys. Res., 85, 8261-8270. 

Phelan, M.J., 1969, Crustal structure in the central Mississippi 
Valley earthquake zone, Ph.D. Dissertation, Washington 
University, St. Louis, Mo. 

Price, J.C., 1977, Thermal inertia mapping: A new view of the 
Earth, J. Geophys. Res., 82, 2582-2590. 


PAGE 30 


Russ, D.P., 1981, The structure of the Reelfoot rift and Its 
relationship to some regional midcontinent features, Geol, Soc. 
Am., Abs. with Programs, 13, 543. 

Sabins, F.F,, 1969, Thermal infrared Imagery and Its application to 
structural mapping In southern California, Geol. Soc. Am. Bull,, 
80, 397-404. 

Simpson, R.W. and R.H. Godson, 1981, Colored gravity anomaly and 
terrain maps of the east central United States, U.S, Geol. Survey 
Open-File Report 81-846, 

Stauder, W., M. Kramer, G. Fischer, S, Schaeffer, and S.T. 
Morrissey, 1977, Seismic characteristics of southeast Missouri as 
indicated by a regional telemetered microearthquake array, Bull. 
Seism. Soc. Am., 66, 1953-1964. 

Stewart, S.W., 1968, Crustal structure in Missouri by seismic 
refraction methods, Bull. Seismol. Soc. Amer., 58, 291-323, 

Strebeck, O.W., 1982, Structure of the Precambrian basement In the 
Ozark plateau as inferred from gravity and remote sensing data. 
Masters Thesis, Washington University, St. Louis, Mo. 

Sweeney, P.H., E.D. Harrison, and M, Bradley, 1977, Geology of the 
Magmont Mine, Viburnum Trend, Southeast Missouri, Economic 
Geology, 72, 365-371. 

VanSchmus, W.R. and M.E. Bickford, 1981, Proterozoic chronology and 
evolution of the midcontinent region. North America, Chapter 11, 
Precambrian Plate Tectonics, A. Kromer (ed.), Elsevier Pub. Co,, 


PAGE 31 


261-.296* 

Watson, K., 1981, Geologic applications of thermal inertia mapping 
from HCMM, 1981 Int. Geosci. and remote Sensing Syinp., IEEE Cat. 
No. 81CM1656-8, 1182. 

Wise, D.U., 1976, Sub-Continental sized fracture systems etched 
into the topography of New England, Proceedings of the First 
International Conf. on the New Basement Tectonics, Utah 
Geological Association, Pub. 5, 416-422. 

Woollard, G.P., and H.R. Joesting, 1964, Bouguer gravity anomaly 
map of the United States, U.S. Geol. Survey. 




PAGE 32 


FIGURE CAPTIONS 

Figure 1 - Digital shaded relief map of topography 5s displayed in 
this image with the simulated sun located at 20" above 
the northeastern horizon. A Mercator projection is used 
for this and the other images in the paper. Numbers 
running vertically on the right hand side are degrees 
north latitude and numbers at the bottom are in degrees 
west longitude. The particular area of concern is 
located in the box bound by 25 to 49" N. lat., and 80 to 
110’ W. long. The resolution of the image, measured as 
the width of a picture element (pixel) is about 7 km. 

At that resolution little structural control of 
topography is evident within the midcontinent. 

Figure 2 - Gray tone image of free air gravity anomalies for the 
same area as displayed in Figure 1. Black areas 
correspond to regions where stations were spaced too far 
apart to allow reasonable interpolation. Bright regions 
correspond to positive anomalies, medium gray to areas 
with anomalies close to zero, and darker areas to 
regions with negative anomalies. The anomaly range goes 
from +66 to -360 milligals. A sketch map illustrating 
structural features evident in the image is show in 
Figure 5. The midcontinent gravity high extends as a 
NE-SW trending linear high, flanked by lows, in the 
upper half of the image. Note the subtle linear low 
extending from a break in the gravity high, toward the 
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southeast, to about 36" N. Tat., 90* W, long. 

Figure 3 - Shaded relief image depicting free air anomalies as if 

they were hills and valleys illuminated by a sun located 

at 15" above the northeastern horizon. Note the 

correlation between the shaded relief of topography and 

this figure in areas of high relief. Note the linear 

gravity low that begins at the break in the midcontinent 

% 

gravity high and strikes southeasterly toward 36" 

N. lat., 90" W. long. A Bouguer image displays a 
similar patten for the linear gravity low. 

Figure 4 - This color image combines a color-coded version of the 
free air image shown in Figure 2 with the shaded relief 
image shown in Figure 3. The effect is to be able to 
discern both the value of an anomaly for a given region 
and an indication of the local anomaly gradient. Red 
corresponds to free air values greater than 20 
milligals; yellow to orange indicate a range from 0 to 
20 milligals; green corresponds to -10 to 0 milligals; 
blue-green denotes -30 to -10 milligals; and blue 
corresponds to values smaller than -30 milligals. This 
figure was the base map used for constructing the 
structural features diagram shown in Figure 5. 

Figure 5 - Sketch map showing major structures delineated in the 
free air images of Figures 2, 3, and 4 and equivalent 
presentations for Bouguer data. The Missouri gravity 
low is a subtle, but pervasive feature that extends from 
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the irn'clcontinent gravity high to the Mississippi 
embayment. The location of the Mississippi Valley 
graben is in part from Kane et al. (1981), A- A' is the 
location of the profiles shown in Figure 6. 

Figure 6 - Topography, free air, and Bouguer anomalies are shown 
for profile A-A' of Figure 5. The profile cuts across 
the Missouri gravity low. 

Figure 7 - Rock types of Bickford et al . (1981) have been coded to 
discrete color values and overlayed onto a shaded relief 
version of Bouguer anomalies. Sun is from the west at 
IB" above the horizon. Color values are as follows: 
Medium blue -• felsic rocks; yellow - granite; orange - 
rhyolite; purple - gabbro; dark blue - basalt; red - 
metasedimentary rocks; green - Sioux Quartzite; and 
light blue - Kev^eenawan sedimentary rocks. 

Figure 8 - Vertical intensity magnetic anomalies with amplitudes 
greater than 600 gammas are shown as black splotches on 
gray-tone and shaded relief images of Bouguer anomalies 
for Missouri. Patterns of small black splotches near 
the southeastern edge of the state and in the upper 
right are spurious and are remnant from the original 
digitization of the Missouri magnetic anomaly map. 
Arrows point to the Bloomfield (lower right) and to the 
Salem magnetic highs (Phelan, 1969; Cordell, 1979). 

Both features are probably due to relatively mafic 
intrusions. Note that many of the magnetic highs are 
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clustered on the edges of the Missouri gravity low. The 
sun is from the west at 15“ above the horizon for the 
shaded relief data. 

Figure 9 - This figure is HCMM daytime thermal infrared image 
covering southern Missouri that has been contrast 
enhanced and transformed to a Mercator projection. Flat 
areas are brighten (warmer) that regions with high 
relief. HCMM frame ID A-A0045-19420-2. 

Figure 10 - The HCMM contrast enhanced daytime thermal infrared 
image has been overlayed onto a shaded relief image 
depicting topography. The sun for the shaded relief 
image is from the west at 20“ above the horizon. The 
overlay tends to enhance subtle topographic features 
while retaining linear thermal anomalies not associated 
with relief. 

Figure 11 - Sketch map showing the Missouri gravity low, mapped 
faults and folds, and 1 inears as mapped from the HCMM 
thermal image, and from the thermal -topography overlay. 
The structural trend is dominated by the trend of the 
Missouri gravity low. Azimuths of structures tend to 
disperse for regions far from the low. 
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HCMM Apparent thermal inertia has been color coded and overlain 
onto a shaded relief nwp depicting topography. Reds corres^Kod 
to liigh A.T.I., blues to low values. 



HCMt-1 daytime thermal infrared brightness is coded as hue, albedo 
as saturation, and nighttime thermal infrared as intensity. Reds 
indicate relatively high values of daytime temperature. 




^ Thi «ad« t* then thnc the ahock vava te> 
aulta tuporteil by Srovm nn4 MacQuiini vV* ifSisirt- 
•il the «>y phaaa ttimaltlgn at Vii CV», favor elia 
choke of t"kon ft Inner >!or« preaaurea. 

tt apt^nra that puru y-tron nlotia aiitlafiea 
the danalty and Vulk ■odulua raquitanenta of the 
Inner core> veochaxUoa erguinentn would favor 
the addition of Ml to pure Irani However If 111 
la edued to y>lroni e lighter aleaent <eueh aa 
nuet eloo be added to eonpenaece for the 
tncreoea In drnalty reauUlnj froi the Hli 
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Shor|i Wove Pet e f or t'yrtta and _lce . Relatin g 

f dJtha ipaTthf arc ora 

THOiWS J, AltRBiS and CllRTSTOPHER R. CREAVEII* 

(Sclanologlcal laboratoryi CaUfornta Inetltuto 

of Tochnolosy, Paeadena, CA 91125) 

(^Preaent nddreeai Aclnntk-iUchfleld Corpi, 
hoe Angalaa, CA) 

Keu ahock and relvaaii wave data In the 75 to 
1850 kbar range for atnglo cryatal pyrltn 
low aptn) danonatrote no deCectlbla phaee cliongO) 
over the proeauco rnnge Inveatlgatedt Combining 
the preeenc d'te with Slnokov cC ali data to 
5200 kbnri flta rjo'iely to linear ahock (U,) 
portlcle (U,) velocity relation UB(kn/aec) • 

5i56 4- LiAA Upi Tlio inforred zero-preaaure 
bulk aound 'pond la cloae to chat obaerved 
ultcaeonlcallyi 5.38 to 5.43 ka/aec and the 
colonae onto conatrnln the Round opeed along the 
Hugonlot CO lie In the range 5.0 to 6.2 km/aec 
froe 384 to 1930 kbar. and conaCrnln the dcnaity 
yarlntlon of Grunetaen ratio to Y • 3.8 (Oo/p)*’^' 
ihla behavior controata markudiy with that of 
Ftto,eS wlilch Initially has In the hJgh-apln 
atace and undergoea an a. 14Z danalty Incroaae 
to a high preeaure phaee below 100 kbar. which 
nay hnvn Fo+* In the lou-apln acnta. When the 
new pyrlto dncn and roconC roonalynea of 
pyrrhofclte or>! u»«d With Bhock vawa data tor 
Iron and autenlc data for the core tho casontlally 
neorly Ideal additivity of the inolnr volume of 
S In Iron eiilfldea can be domonatcnCed. Depend- 
ing on thermal protlloe aaaumed. aono 9 to 13 and 
10 to 13Z S nro Inferred In tho liquid outer 
cure from the pyrlto and pyrrhotlco data, 
reapectlvoly. 
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nvna nlC B nf chemically En rtchod, Int erp ally 
Heated Plumea In the Hantle 

Henry R. Eby and David A. Vuan (Dept, of Geology, 
Atleona State Unlvaraity, Tcope, AZ 85287) 

It haa been propoaed chat plumea are ftnleo- 
anplltude nanifoatattono of locullacd thermal 
Inatabllltlea produced at the lower (hot) thermal 
boundary layer, owing to the deatablllzlng 
effects of Btrongly comperature-dependenc mantle 
rheology. If the prlnUlve deep manclo lo en- 
riched In tndlenctlve heat toutcca from Inhomo- 
geneoua plonototy accretion (Andarson. 1975), 
then thermal Inatabllltlea moy fumlab a viably 
meebanUm to axtract the enriched hoot sourena 
Into tho npuelllng, which la then propelled 
rapidly upwards by thcee moving heat sources. 
Accordingly wo have modelled tlie thermal- 
mechanical atructutao of thaae chemically en- 
riched. subaolldus plumea by uenna of boundary 
layer techniques In which both Intornal heating 
with a locallaed apatlal variation nnd a strong- 
ly temperature-dependent ylacoalty, eluirnctnr- 
Isclc of mantle materials, ate Incorporntod. 

Ii;e equntlona govamlng the conservation of mass, 
energy end moraentum nre solved numerically* Two 
different modes of heating have been studied* 

(a) the plume la dtlvan aolely by a conccntrncod 
amount of radloganic heating (b) the energy 
source la derived from both the localUed radio- 
activities nnd a hot adlahnr.lc tempetneuro Im- 
posed at the ccntBrlino of the plume. Wo find 
that for mechanism (a) to produce velocities 
greater than 10 cm/yr requires radioactive 
concentrations equivalent to those typical of 
granitic rocks. On the other hand, narrow, 
fnnt mantle plumes, 0 (n/yr), are easily achlovn- 
hle tor heat production rates clinrnctorlstlc of 
oceanic basalts and adiabatic temporature eg- 
ceases of a few htmdredn of degrees. 

Anderson, D.L., Chemical Plumea In the Mantle, 

Cool. g oo. America gull .. 88 , 1593-1600, 1975. 
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TlmfDependent li'W.srkal Hodaln nt si ngle and 
iio^'l erj g/»r rMan t'i ejCiinv Vrtln n 

A*^JL* <">8 *• (Dapt, of Terrestrial 

NagneTlan, CarnaqU Institution of Wash (niton. 
5241 Proad Branch Poad tlW. Wathlngton t>G 20015) 

Contrary to veochtnlCAl evldenca. current models 
at puintle convection aceumlni stendy state a<otlon 
favor alnile-layer convection. Rslaxlni th« ataady 
state aaeuntition, we have studied the tine evo- 
lution of two-dlneiialonal Housetneaq convection. 

A survey Involvlni variation of aoveral paramv* 
tera has been completed. Each model hne Cartesian 
qaometry. exponential vlacoalty dependence on 
temperature, nnd frea-ellp boundaries. Initial 
conditions Were chosen as a fluid at rest with a 
uniform tamperature gradient (pure conduction) . 
randomly perturbed at the Si level. Paronatore 
varied Inoludei choice of Initial noise. Ra 1106 
to 108),aspact ratio A (1 to 2), heatlni source 
(all bottom, all Internal, or half and half), and 
s|Hitlal variation in the temperature along the 
bottom boundary (2St amplitude cosine variation). 

The time development if ennaitiva to the choice 
of the Initial nolae, except when varlatlcns In 
the bottom temperature ere present, which then 
determine tha tranalCAt call structure (pltunai. 
riea over oach tempereturo maxtmien). For Ra«10', 
the Initial ctnile-layer ntructure focna a tran- 
sient phase of double-layer convection, with as 
many as 4 cells In each layer. *teady state lolu- 
tiona In all cases favor sinile-layar convection. 
Inorsaalni A Incrsnaes tho tendency foe the 
double-lnyar tronaiant to occur. Varying tlie loc- 
ation of tha heat source has little qualitative 
etfeot upon the evolution. Double-layer etructura 
occurs without the targe temperature difference 
between tha two Isysra required in steady atate. 

The timescale for existence of tlie double-layer 
strueturo la 1 b.y, (with Ra>107, os appropriate 
for whole-mantle convection). When compared to 
the timescale for plate motions (about G.l b.y.), 
t)ie results imply that tinu-dependent modela are 
neceiaary for understanding mntlo convectlonf 
transient periods of doublo-loyor convection may 
have oxistad In tha past, or may exist now. 
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T21D-1 

Gravltaclonally Induced Strosses In a Thin, 
Spherical Shall 

A. K. COODACHE (Division of Gravity, Oeothermlcs 
and CcodyunmlcB, Earth Phyalcs Branch, Ottawa, 
Ontario, KIA 0Y3) 

In situ stress neasuremanta generally Indicate 
that the vertical eompreaslve stress la aero at 
the Burfaee of the Earth and Increases with depth 
approxlniately as the overburden pressure. On the 
other hand, the mnan horlsontal comprcsalve stress 
la non-xero (nlO MPa) at tlie surface and Increases 
with depth less quickly chan the overburden prea- 
aure does. The rnault is cliat the ratio of 
horirontal to vertical across la geeator then 
unity near the surface ona only tends to foil 
below unity at dapths greater tiinn 1 to 2 km. This 
bebaviour (in particular the presence of non-zero 
horizontal compreastve aCressea at Che surface) is 
prcdleced for gravicacionally induced BCreasea in 
a thin, homogeneous, clnscic spherical shell en- 
closing n ''fluid" under n prasaura equal to cha 
total weight per unit area of the shell. The 
horizontal compreastve etreas at the surface la 
related to the thickness of the shall. Assuming 
that only a rhln portion of the Earth' a outer 
cruat can sustain shear neresaes over long periods 
of geologic tine, data from Canada and southern 
Africa indicate that the oflectlva thickness of 
die clnscic crust In theae nroas la of the order 
of 4 and 2 ka respectively. AUhnugh two-dlmen- 
sloR.*!! elastic theory does not predict finite 
horizontal compressive scrosses at tlin Earth's 
eurfoco, die results of the thin shall theory can 
be used to help correct tho calculations nf 
gravitationally Induced atressea In simple two- 
dimensional geological models. 
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The variation In heat (low in Arlionn lias btsn 
Intatpreted utilizing n variety of geophysical 
end gsothefitcAl dsin, The depth- to-Curie Iso- 
ihem has bien previously determined in Arkons 
by SF dtral ».*iilyils, spectral M'chlng, nnd by 
exaslntng the general amplttude varktions of 
aecomainatic anomalies, there are aig.alfkant 
dtBsgrow.Mnte between the interpratationa of 
much of dies* data. In th,s study, the Curie 
depth haa been annlyzad by inversion of ovar 60 
nagiietU anoaalUi, with the gackua-CilbeCt 
technique providing etatlstlcal anstyala of the 
raauUa, Vertical rasoluelon Is +2 kn at 20 km. 
Depths range from 5 ks below asa level to 21 km 
b.a.l. tho reaulta agrea In RSnetal with spac- 
tral matching ratulci, but not well with the 
otiicr typae of Curie dupth analytes. Sliallow 
depths ate Indicated along the trsnaUlon bo- 
tween the Colorado Plateau and Batin and Range 
provtncei, In the White Hte., In aouthwest 
Arizona, end in zonthansc (Tucson area) Arizona. 
Central and aouth-cencrnl Arliona In the Basin 
and Range show anousioualy deep (20 k.m) Curio 
depths. An nssumpeton la made that the thick- 
ness of magnetic bodies la deternlnad ctiensally, 
not lithologically. Although the lateral varta- 
tlon In heat flow oeasurements genarally agree 
uUli liitersl changes In Curio depths, few eois- 
pletoly corrected values are avatlnblo. Some 
correlations exist with crustal seismology ex- 
periments. FraUmlnary Hf ciodeUng by K, A. 
Ander of Los Alamos Kaclonsl Laboratory Indlcatm 
extremely close correlations with ratactve 
variations of Curio depciis In various parts of 
Arizona, Kumorlcsl heat flow models using the 
Curie temperature as a boundary condition have 
been constructed to Conpare with hent flow and 
KT Interpretations. 
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Linear Basement Structure In the Hldcontlnent, 
U.5.A. - Gravity, topography, remote sentl no, and 
seismicity data , 

RjE. ARVIDSON. E.A. GUINNESS. J.H. STfiEBECK 
{RcBonnelT'Center for thp Space Sciences, 
Washington University. St. Louis. Mo,, 63130) 

Soire 500,000 Individual gravity readings for 
the continental U.S. were reduced to free a1r anJ 
Bouguer anomalies. The anomalies were then fil- 
tered to produce maps with 255 contour Intervals 
using the algorithm that generated the Pioneer 
Venus altimetry maps. A linear gravity low, 
with an amplitude of about 30 mg. extends from a 
break In the midcontlnent gravity high In S.E. 
Nebraska to the Pascola Arch (horst) In the 
Mississippi f.ibayment. There Is a suggestion of 
an extension o the northwest Into Wyoming based 
on 0) The alignment of a sharp gradient (but 
only a few mg A) In the free air data along the 
North Platte River, and (2) The alignment of 
fractures In the Prccambrian of Wyoming. In 
southern Missouri, as shown by digital topography 
and Heat Capacity Mapping Mission data, the 
gravity low seems to nave controlled a number of 
basement normal faults, In addition, most of the 
New Madrid seismic epicenters occur In the crustal 
block defined by the Intersection of the low with 
the Mississippi Valley graben. The origin of this 
linear basement structure Is problc.-natlcal, 
although the fact that tho structure cuts across 
age and province boundaries Is suggestive Of a 
transcurrent fault. 
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Linear Structures In Missouri Seen with Heat 
Capacity Mapping Mission (HCHH) Data 

J.W. 5TREBECK . E.A. GUINNESS, R.E. ARVIDSON 
-^RcBonrioTT'Center for the Space Sciences, 
Washington University, St. Louis, Mo., 63130) 

An apparent thermal Inertia rap was produced 
for tho state of Missouri using HCMH night/day 
thermal Infrared (IR) Images and an albedo Image. 
Differences In rock types (between Prccambrian 
Igneous rocks, Cambrian section, and younger 
Units) and some structural linear features pan 
be seen on the map. These structures were empha- 
sized by overlaying the thermal Inertia images on 
digital shaded relief maps depleting topography. 

In general, there Is a good correspondence between 
structures seen In the HCM.M data and mapped faults. 
Often, fICMM structures extend beyond limits of 
mapped features. In addition, a number of natural 
linear features were detected that do not corre- 
spond with napped structures. Finally, tho ob- 
served trend of many of the HCWI structures Is 
coincident with a 30 mgal, 140 kn wide, Bouguer 
gravity low that traverses Missouri In a tW-SE 
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McUoniirll Oiiier lor ihr Spjir Stirnccs 
Wa^hiiiKion L'mvei«ii> 

Si l.ouii, MikMiuri 

Introduction 

A (liKital iu|Migr.iplii( iiiloi- 

niaiion and j{ra«i(> data cxIm that rmrr llir 
conliiiriital I'nilrd Slain Koi iiiMaiiic. dal i 
irii dial can be- pnrrhaxrd Irom ihr Naiional 
(kranir and Aiinotpliriit Adininixlialion 
(NC)AA) includr (n) cirvaiiont axrraxcd mci 
^0 i ol laliludr and liingiliidc for ihr coirr- 
niinuk I'niird Slam and (b) uvrr a hall mil- 
lion Kiaviiy iradiii((x rr-iin iiaiionx dixirihiiird 
aciox% ihr rniiir country, alihoiiKli in an ir- 
icKiilar inannrr. Hoih ihr io|io)(raph> and 
ihr )(iaviiy data xrit air lundamrnlal lo ihr 
ui .I^-ixtandiiiK ol ilic xiiiuiiirr ul ilir ciuM 
and lilho'phrir, and rx|>rcially lo llir under- 
ttandiiig ol llir irlalionxhi|ix l>rlurrn Io|mi)(' 
raphy ami \liiuliiir |.Wrf>inwu el at , l'J7'J: 
MfXult, 1981, and oihrr»|. 

I hr inirnt of lliix woik ix lo chow ilir ap- 
plkahiliiy ol lairly Mandard imagr-procrcciiiK 
Irthniqiirc to piiKrcciii){ ,ind analyziiiK large 
geologic data sets. In pariiiulai we hase uii- 
ii/rd iiiiagr-lilirring let liiii<|iiec to iniripolair 
lirlwreii gravity clalioii liKalionv in older lo 
piuduce a rrgiilaily tpaird data array that 
prrvrrvi-v driail in airav vvilli giMKi coverage 
and dial piiKliicrv a roiiiinuoiiv tone image 
radic-r ihan a ronioiir map W'r have used 
viandaid iinagr-piiKcvving lei lini(|iirv to digi- 
tallv irgivtrr and overlay lo|Migraplii< and 
giavily d.ila, and ue have displaved die data 
III wave that rnip!iavi/e viihilr hiil |>rrvacivr 
Mriicliiral lealiirev In diiv pa|>rr we divciivv 
die lechiiiipirv lived, hrieHy deviiilie dir 
piiHiuciv, and illuvtrair dir |Mitrmi.il ol die 
medicKlv hy divciiwing iiihilr linear striiciuirc 
dial ap(>rai in ihe picKrvved data liriwrrii 
die midi ontmrni giavily high and die Appa- 
lai liianv. 

Data Processing 
Techniques 

h.ila piiKevvmg uav condiicled on a l‘l>P- 
II .M mmuompiilrr viidi inlerailivc image 


OF POOR QUAUPt 


Moiage and divpiav |>riiplirralv Soliwair iv 
havrd on mmi-VI(^\!<, widi a vaiiely ol ap- 
pliiaiioiiv programv locondiui lilirring. en- 
haniemrnl, and gc-onieirn operationv I hr 
reader iv lelerrrd lo viich viaiidaid irlri- 
enerv av .W»ik I l‘i8U| lot liiidiri inloimaiion 
alioui general trilinii|uev lived in miagr pim - 
evving Our viliwair wav drvrlo^red patily ai 
dir |et Piopiilvioii l.alKiiatoty; p,ird\ ai die 
I' S. (•rological Survev. Magvtali, Arizona; 
and pailly hy iiv. 

For piirpovev cil ihiv pa^iri die dynamic 
range ol ImhIi die io|M>giaphy and die gravity 
data vvrrr htvt vcalecl to hi widiin the tange 
occupied hv a hvir vaiiahle (i r., 8 hiiv or 2.’ih 
divciele valuev) I he hy ie-eii« mled d.ita vveir 
then vioied in d..ia arravv, where dir aiiay 
riemrni tpacing wav SO v m holh Lililudr .iiid 
longiliidr For the lofxigiapliv ihiv picHediire 
priNluced a regular array of hyie rnciKled 
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dala dial had 28HO lowv llaliliide) ..nd *080 
coluimiv (longiliidel I hr dala airavv covrieil 
a laliludr range horn ?.'•* lo 19 ,\ and a loii- 
glliidr laiige lioin l'8i* lo I'.’a'W. 

Ihr levidling lo|Migiapliv dala aiiav could 
hr divpiavrd av an image ll die hvir-riucKlrd 
d.ila weie convened lo color or hrighinevv 
valurv An alleinalivr and i|uile iiveliil rnrlh- 
inI ol depicOng lo|Migiaphv it in generate a 
digilal, vliadrd lelicl map figure I vhirwt 
viicli a producl. winch veat grneialed on a 
vrrvioii ol die dala lilr dial wav liitl liaiiv- 
loitiic'il lo a Mercaicn piojei non, digiiallv le- 
duced III vi/r hv a l,n lor ol 7, and dien lived 
lo gener.ile a reliel map I he vhading priMrvv 
iv haved on die avviimpOon ol a pholonirliic 
lull! lion loi die VIII lace (Fomnirl Seehger 
luriclioii III out cave) I hr hiightiirw av- 
tignrd lo a given aiiav rlemenl IrHalioii dr- 
prndv on die local vlo^ir magniiiide and di- 
irction and it lelaiivr lo an avvunieil volar az- 
iiniiih and rirv.ioon [Hatuni rt at , I ')7j>|. 

Ihr gravity dala were reduced lo free air 
and lo Hougtirr aiionialiev cuuiiriv of d.c 
Gro|>ovitional Drpaiinirnl. Peleiivr M.ippmg 
Ariovpacc ,\grncy Onirr, St l.oiiiv, Mivvciii- 
ri I hr irrrrriite held at vea level dial wav 
uvrd III die rrduclion wav haved on die I'.8i7 
Inirrnaiional (iiavitv Formula, with all die 
dala ref erenird to die 1971 Inirrnational 
(•lavity Siandaidizaiioii Net Free air aiioina- 
lirv were compuird liom the following for - 
niulation: Sg, v • g S08hh — >, v*hrre 
•^Jifi V ** fire air anomalv (nifial). g • rr.nl- 
ing, h rirvalion (|K>viiive down to groid), 
and > w iheorelical giaviiv value. A vccond- 
oider tcim v>av added lo the elevation coriec 
lion when die magmiiide ol the coiieclion 
exceeded 9.1 m(•:d Bougiier anonialirv wc ir 
compuird on ihr haviv ol a vl.ih with drnvilv 
‘i.l>7 gni'em', living Ajf, ^ - 9 1 1 I9li, 

whrtr Agiiv ^ Uougiiri anomalv in m(•.d 
L.ck,iI ierr,iin coireclioiiv were not included m 
die Bougiirt .nionulv i onipiiiatiiiiiv. 

(ienrr.iling imagev lioni the gravitv dala 
pirvrnlv a mote cotnplicaled piuhirm dun 
die lopogr.iphir dala vincr ilir giaviiv viatioiiv 
are not l<K,iird along a legiil.ir grid For ImiiIi 
die f ler air and Bougiiri aiionulirt die d.ila 
wrrr hivl vcaled lo hi vtidiin a hvir variable 
(i.e., 2.'>li divcreie value-,). Ihr vcalecl dala 



Fig f Shaded lelirl m.ip depicimg \().\ \ .S9-v elevation aver.igrv lol llir ronihient.il 
I'nited Statev. .Niimlirrv running veiticallv along dir videv ol the ini.ige air degreev ol 
noidi laliludr, while tiumlieiv along die Imiiioiii .iie degierv ol vtevl longiliide I hiv .md 
other mapv m the p.i|>rr arc Mercaloi pio|eilionv I he vnnulaled vun wav :a*i at .91’ alHive 
the wevlrro hoiizon It.itk Mch kv aie.iv aie niiwnig dala 
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f'tg. 2 . Sludril rrlirl iiup ilrpkliiig lilirrrci hre iiii llljik .m-j« iiulM.ur rr- 

gK)ii« coiiuiniiiK riihrr >m (1 (Iju ur kmi liitlr tijij lo luvr vuiiuit jity \ jIkI iiiirrptiljiioii*. 

I lir kimuUtcd »uii kci at IV aUne llir Mrurrn liori/oti VaiLiiioiik in ihr anoni.ilio 


arc. Ill cflctl, illuiiiiiiaied at d they ucir a tcl of loitograpliir ndget and talirtt. for exam- 
ple. the niidioniinenl graviiy liigh (-IV.N’. ‘.»r\V). uliuli it a |iotiinr anomaly. Maiidt at a 
ridge, vthile die gieai talirt ol (aldornia (37*.\. IVO“\V). ttlmli it a negalitr anomaly. .i|h 
prait at a iioiigh 
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tierr plated into an airay ttidi die tame rlr- 
iiirni ipating (3U t) and grogiaphit (otcrage 
at die lo(Mtgrapliy ilala air.iy. I lie rirmeni 
liMaiHm (loteti lO a gitrn gratiiv tialion Itaa- 
iHin Hat attignrd dir ttalion taliic. ilieirliy 
pKNliMing an airat rxtiipicd in |>ari liy real 
data and in pan liy blank lonet 

\Vr ihote lo ute a loiurpliiall) tiinple. IhiI 
Mivtrrfiil. filiering approarli lor inter |Milating 
leiiteen ihr leal giatily data (xMiilt I he 
irilini<|iic wat drtelo|ird lor piocrtting lunar 
(oiitortiiim data •.ml ttat alto iited lo pio- 
duic aliimriiy niapt fioiii dir appioximaiely 
IIHI.IMMI rirtalioii rtiimairt ulMainrd lit die 
I’KHirei Vriuit radar mapper (f./ioJiM and Sn- 
dtrhiiim, 11)77. /V«r«gi/f r/ <i/.. I(IM0| Ihr 
Irthniqiir imoltrt iitc ol a tp.iiial filtei dial 
genrrairt a moting atrragr lliiongli ihr data 
array. Ai any gnrn potiion ttidiin die data 
array ihr filler oetiipiet a Inix of .V x .V rlr- 
nirnit I hr inraii tahic rd die giatily aiionia- 
liet for kialiont tiMaleil ttidiin ihr Imix it uteri 
at ihr inierpolaied tahir for the box irnler 
only if an anomaly value doct not aliraily ex- 
it! at dial liMalion Shilling the filler oter .hr 
airay aiirl ir|>raling ihr o|>rralion foi caili 
eirnirni ItKalion irtiilit in a pai dally intrrpo- 
l.iird data tel. blank airat tiill irmain in 
ruiirt ttiihoni any dala anil in /onet ttidi m> 
liillr dala dial a iiser-definrd iialitiital 
ihirtholil wat iioi inel I he ihrrthold it ilr- 
fined III lernit ol the lr.i<iion ol eirmrnit 
ttidiin the Ixix dial niiiti hate t.iliil data in 
order lo tepl.ne die mid|Miint eleiiielil with 
ilir ateiage. .-ittiiming dial die midpoini did 
noi aheadt hate tahd rlaia. 

In pi alike die giatilv dala were piiMewiI 
hi itting teter.il Idler widdit. iK-ginning with 
a 3 X 5 rirnirni filter and ending with a *il 


X 21 eirmrnt filler. I hr ilnrthold wat tel 
tiirh lhal al leatl 20'a id ihr elrincnt hwa- 
liont at any gitrn filler |M>tiiion nniti hate 
hail real data tahirt lo lie iniri|M>lalril I he 
choir r ol filler ti/rt wat gotrinril by die M.i- 
lion tpatingt. wtiiih atri.igerl tetei.il kilome- 
Irrt but tariril fioin Inindierlt ol ineiert lo 
teteial lent of kilonielert. ,\ .3 x .3 elemeni 
filler, for intlanic. cotrri INI t of an in ImiiIi 
lalttiidr and longiiiidr. oi alioni 3 km along a 


line ol laiiitidr. .\olr dial legiotit with tela- 
lively tmall ditiaiuet lieiwren giatiit taliirt 
were inirr|xi|.iirrl unit oter irlatitrit thoii 
ditiaiurt, tiiirr once an rirmrni wat attigned 
a rlata t.ihir ii w.it noi .dinirri bt laier Idler 
pattet. f hr net etlrri of the filiering o|>rra- 
iHint w.it lo geneiair a irgular aii.it of giati- 
ly rlaia. wilh original riala lift iniart and wiili 
iniri (Mil.iied il.iia iN-iwern die oiigin.il t.ihiet. 
.Some areat in die ariay. even wiih a 21 x 21 
rleinenl liltri. tiill h.id Iiki hide dai.i to allow 
miei |Milaiion. I hote areat apfieai bl.irk in all 
of oiir ilitplayt. 

'I he liliererl giatily anomaliet could l>r dit- 
pl.iyed at giay lone oi at color-ciMlrrl iinaget. 
At with loyxigraphy. we hurl dial a very 
rlleri.tr titual ditplay it a th.ided trlief map 
Figiiirt 2 ami 3 rlepiii tfi.nlrrl relief ma|>t foi 
die f ree air ami Hoiigiiei anomaliet. retyK'i - 
litely. (aimpatiton ol ihrtr pnMlnclt with 
pirtioiitly ptihhtlird contoui inapt (r.g.. 
Wnolard and Jnrsting, I'.Ni'l; .WiffOoio fl al., 
Il)7‘l| thow timil.ir hioarl-tcalr p.itiernt 
llowetrr. die inapt thown in Kignirt 2 .mil 3 
ditplat inhrrrndt more inlormaiion. The 
rra tout are ihirrfold firtt, a linn giid t|>,it- 
ing w.it ntrd to iniripol.iie iK-iween ttalion 
readiiigt ih.in ant piihlitheil map wr hate 
been able lo exaininr. Second, dine are at 
many tontoni inint.ilt at iheie air tablet in 
a byte taiiaMr (i e . 2.M>) Most lonioiir inapt 
lypirally hate a do/ni or to tonioni inieit.ilt 
lliird. arr.it with .i l.iige inimijn ol ilotely 
ip.iieci tiaiiont iri.iin det.nit id ihr anomaly 
p.ilinnt 

'Ihr Idiriril giatily anoni.ihrt ran alto Ik- 
chgilallt oteil.itrd onto die to|Migiaphy loi a 
titual ditpl.it ol the I Ol lel.iliont Itnwern die 
Iwu taiialilet. I he In hniipie dial we used 
wat lo liitt conteri die htte tali r loi r.iih el- 
riiirnt in the gi.it iiy air.itt to a given loloi 
ifiguir-ll The I oloit were c hotel intmh.i 
niannrr dial blue toriet|M)iided to ii<e moti 
nrg.ilite anoin.diet .iiid led lo ih'- iiios. yxiti 
lite anoin.dirt Oi.inget. tellowt. am! gii-rnt 
weiechottn lo leptc-teiil mieimeihale v.diiet 
( ailoi -tiHled imaget dipicling the anom.i.t 
palleiiit weie dtgil.dit oteilated onio a sh.id- 
rd leliel map of the lo|Migiapht timpit In le- 
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hg. I. Shjclrd rrttrf itMp ilcpttiiiiK filierrd bouxucr aiionuliet 'I he map wa« pri>- 
duteti in the uinc manner at ilic Free air anomaly map in Figure 2. No l<Mal terrain cur- 
reciiunt were incluiled in computing Buuguer anntnoliei. * 
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hg. 4 I hit diagram thowt the appiuxi- 
male coloi -ending Mhenir utrd to produce 
the cohir contmiti fui the gratiit aiiomalt 
tnapt thovin on the cover ul tint ivtuc o( £o>. 
Hittogtaiiiv thowing the Iraciional atea ckcu- 
pied h) euh anomaly value tiere uved to at- 
sign cnlort in tuch a manner tiiat each color 
value occupiev an area t otnparahle to any 
other color value. Alxnit -l&divcrete color val- 
uet (i.e., ditc retr hue and latiiration valuev) 
vvcre iivetl in each anomaly map On the oth- 
er hand, the vhadrd relief map vertioiit of 
the anoinaliet each toiitaiti 2jh vhadet of 
Rf»y- 


placing the hiightnetv componetit for each 
color with the hrightnevt from the torie- 
ipondiiig IcMation in an art ay containing the 
tliadrd lelief map. I he tiet retult it a color 
map where the hrightiiett it nicNiulatrd by 
the th.icled telief map, while the color hue 
(domitiaiU wavelength) and tatuiation (de- 
gree of purity) are controlled by the giavity 
anomaly valuev. Note that any color can be 
uniquely defined by a hue (dominant wave- 
length), tatiiratioii (puiity), and biighttiew 
\Slatrr, Ui73) Such picKl'JCtv ate thown for 
both the filtered Free air and liniigurr anom- 


aliet on the cover of thit ivtiie of £u In both 
catet the inapt were tranviormed to a l-ferca- 
tor projection, digitally reduceit by a foe or cd 
7, atid then prcKetted at ditciitted af>ove. 

Discussion 

A tketch map thowing tome of the tnajor 
ttructuiei evident in the color iiid vhaded re- 
lief niapv it thown in Figiiie 5. Foi thiv paper 
we will rettrict our diMiittion of the pnH- 
etted data to a pottion of the niidcontinent. 


where we have a continuing interevt in the 
cruttal ttructure and the telationthipv be- 
tween batement frocturet and ore depoiki. 
Such a dinuttion alto tervet to illuttr.xte the 
kmdt cd information that can be gleaned 
from tofiography and gravity data that liavt 
been priHetved in image format. 

I’he d.ithed linet in Figure S define two 
major tttucturet in the niidcontinent: a newly 
ditcovrred feature we c.ill the Mitvinri Riavi- 
t) low and the Mivtivtippi X’allev gialwn, av 
defined by Kunr ft itl (UIKI) on the hativ of 
gravity and acroniagnetK aiKimaliet. 1 he 
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AlitwMin Ki4\iiy low 1411 In* w-rn ihi ilir («4or 
<ind khiMlnl irlirl iimim jt » I lO-km-wKir Irj- 
liirr ihji hr)(in% Ml m IhcmI nr.il llir Miiiilirtii 
rilifr ol ihr nii«Koniiitriii K*>*-i*y 
r\trii<U m/kim Mi««oiiii Mini iiiio ihr Mitti* 
M|>|H rmiMinirni In ilir ( 11111 . 1)1111111 ilir low 
lMkr« llir loini ol m Iioi%i IiIin k (.illnl llir I'm« 
(oIm Aiili |/'.n’i*i ami/ A/rfiinnii, I07&) ‘ihr 
low ha« a IVougiirr Min|ililiiHr ol mIhhii 30 
ni(«al foi llir legion jii%l 10 (hr noiihwr«l ul 
(hr riiilM) iiU'iil I hr niMgiiiliiilr of (hr jiioni- 
aU (UKKrtIt M (Irrp Uiw iiiriK MiuKliir III 
fan (hr aiioiiial) i« (onti»iriil wilh dii iii- 
(rra«i (iii»(al lhukiir«r ol 3 km al (hr Molio 
or wilh a cruti dial hat a ilriitilt of 0 I gm/ 
cm* Irtt (hail llir tiirioiiii(ling airat lor (hr 
lirti 4 to M km Iwlow (hr tiirfatr |.S/o4«A, 
I9H2| Analvtrt of iriiioir trntiiig data and 
ground ktudi(*t loi airai o\rrl)iiig (hr giaiity 
low drnioiitlralc dial (hr irgioii.ii fiariiiir 
irrndtair (oiiHidrni wilh ihr tlrikr of lint 
drrp liatrmriil fraiuic, tuggrtimg a fail dr- 
grrr of innlrol ovri llir |iaiirrn ol fauliiiig in 
(hr aira [GuinHru rt ul , l'JH2|. 

1 hr iiitrricciioii ol (hr Mitwiuii gravity low 
and ihr MitMt»ip|n Vallr> gr.ilirii hat drai ly 
trivrd to lonirol ihr rm|il.i(cinriii ol two 
malic pliiloiit drfmrd by IhltUuhratui ft til 
(1977). (hr ItliMiiiiiiigioii iiiiruiioii 10 dir 
noidi anil dir Caiviiigtoii iiilriition to ihr 
Miuth I l.r iwo iilnioiit can lir vrrn at |Kiti- 
livr aiion*alirt in dir Irrr air coloi map and 
at bump* on dir thadrd irlirf ni.i|» Adili- 
lionallv, hgnrr (i it a gray lour iniagr of dir 
hlirieil I irr aii daia that hat iM-rii lonliati 
ciihaiicrd to thow (hr vaiialioiit in die mid- 
cotiiinrni. In tiiili a display dir nioti nrgaiitr 
anonialirt aic daik. and lett iirgaiivr atioina- 
lirt arr biigln 'Ihr Mitwmri giaviiy low it 
(Irailv diHrrnibIr, and dirrr it a tnggrtlioii 
of a N'h-S\V (iriiding tmicinrr dial, in iail, 
oullilirt dir Mitti**i|ipi \'allr> gialirn. llir 
wliiir ilutt irpirtrni mrr lOltO t'ardi<|iiaki' 
rpiienirtt irroidrd bv dir Si I iitiiv I'nivrrvi- 
IV i(.‘itinii nriwork iM-iwerii 1974 and 1979 
|.Sr<iui/rT ft al , 1977). CIrailv, dir cpiirinriv. 
whkh driinratr dir New Madrid vritmiiiiv 
high, arr loiiiriiliaird in d:r iruvi.il blm k 
drfinrd by dir iiiirrmtioii ol dir Mittniiii 
giaviiy low and (he Mittitvippi \'allrv gialirn 
Thiit ibr inirrtrtiioii of ihr«r two crtivt.il 
tliuclurrt trrmt lu have piovidrd Iwidi a 
zone of wrakiirtt for dir rmplairnirnt ol 
plutoiiv and a focut lor irlratc ol ttrain riirr- 
gv via raiihqiiakrt. 

I \ammalion ol dir lire air color map tng- 
grtit an rxirntion of die Mittonri gravity low 
fioni dir midconiiiirm giaviiy high lo dir 
iioidiwrtl, |K‘ihapv at fai at (o dir big Horn 
nph!t III Wyoming. I hr tnggrtiion it hatrd 
on an aligiimriil of a thaip birak in the Krrr 
ail anmnalirt with the .North I’lattr Kivrt v.il- 
Irv ill Nrbintk.i and with fiatniirt in Wyo- 
ming 'I hr iralily of vnch an extriivioii iiiiivt 
irinain o|>rii to ipirttion until riiilhrr tindy. 
II the rMdition it grnriiially irlaird in ihr 
Mivwinii graviiv low. the lombinrd Iriiglliv of 
ihr IWO Iralnirt would br ovri l.'itMl km 
'riirir air other ttiiuuiirt of conipaiabir 
length on ibr r.nih. iianirtv ti.nitctnieiil 
faulit tnili at dir Alivn I agb of Cdiina |.W»I- 
rjir mill I'a/i/mnifr, 197.5) It i* loiuriv.iblr dial 
the .5|ivtonti giavitv lot* and dir rxirntion 
into Wyoming are pail of a liantr iirrrnt lauli 
tvtirm. tetiiontof tthii h havr brrn iractivai- 
rd dining vaiioiit tiiiir iirri'Mlt Siiih an ori’ 
gin It atwi (oiiMttrnt wilh llir obvrrvaiion ihat 
ihr til lit loir riitv .most the batrnirni .igr 
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fig 5 Skruh map tliownig tomr of ihr majoi tliiidnirv tern in llir giaviiy aiiomalv 
inapt pirtrnird in lint p.qiei I hr Mitvinii giatiiy low can lir tern Iwti on dir tha>lrd ir 
lirl map trrtiont of llir aiioinalirt I hr Iwo pluionv .nr (hr Idiximfirld pliiion lo dir north 
and ihr Covingion pinion 10 dir tmiih I brt can lir trrn at iril dolt on ihr Krrr air color 
ma|>t and a* biinipt in dir tliadnl rrlirl anomaly inapt. I hr liliMimfirld inirntion it liwaird 
al ahoiil 37*N, 90*W. Tlir Mittittippi tallrv giahrn liHalion it lioni kanf ft ul. |l9M|j, al- 
though (hr oinlinr can lie tern in ihr thadrd irlirf aniMiialy m.i|»t. 



hg 6. I.ray lour vrrtion ol the lihrird firr air anoiiulirt I br map Ixiimdaiirt rxirnd 
fiom 34’ lo 43*N. 87* (o Wi’W I hr iw > hoii/ontal linrt conrtpmd 10 37* and 40*N lali- 
ludr. while ihr two vrriicnl linrt air 91* and 9.5*W longtoidr In diit map. biigbi airat b.ivr 
large I irr .lii anoin.ilirt. while daik arrat havr tm.ill (i.r.. highly nrgaiivr) valiirt. I hr 
white dolt air raMlwpiakr rpurntriv lor llir prrwxl fioni 1974-1979. at motdrd bv dir Si 
I oiiit I nivrrviiv tritmic iiriwoik. Ihr bliMimfirlil and (!ovii<gion pliiionv arc dir biighi 
ciiciil.it frainrrt |iiti lo die iiorih and toulh. rrt|>r« iivriv . of ihr rpitrniri ihitirr. I hr Mit- 
touri giavity low rxirndt lioni a break m die mKlionliornt gravity high, ihioiigli Mittonri. 
and ai rott ihr Mittittippi vallrv gralirn. Ihr iniritrtiiont ol dir giavitv low and (hr gia- 
lien air die hxaiiont of dir niliiitiont. In addition, moti ol dir trivinxiiv it « mu rnli.itrd 
within (hr crioial Idixk driinrd bv die micrvciuin ol die giaviiy low .md dir gr.ilM-n 
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•imI |/i»vimr iHHiiutjiir* ilrfinrd by tinkfutfl 
rtm! |IOHI| in MitMMiri ^or inMunir, ilir Ir^- 
iiiir loiiUl luir loiiiinl duiin|( ilir |m>miiI.ii> 
ni (ollitioii %«riil ibjl (iMHliNril ilir (•irnvillr 
uio|(rny \Unrry ami lluikr, IU73| 

Implications for the 
Earth Sciences 

‘[hr ruiiiir will »rr aiti iiuiriiird 4«ailjlHliiy 
Ilf * minilirr u( iiiliri Uigr, ^niKtjuliMailly 
oiirntrd d^U »rii ilul arr applKaltlr In rarih 
Miriiir piiililrini. [oi iiiMaiMr, ai pail of ihr 
NalMMul L'rjiiiuiii Rrtoiitir l.valiiaiion l‘io- 
Kraiii, ihr Drpariniriil of Liicrgy hat at- 
quirrfl aiihortir diKital Kamiiia ray (iiri laifMi- 
ailivity) dau, arioniaKiiclii anonialir*, iiiiil 
^ydro^r 1 HllrlllMal data for moti of ilir loiiii* 
nrnial I'liiird Siairt. Siuh laigr data »rtt air 
i|Hly tuiicd for u«r with ihr kind* of irrh- 
iiii|iir» drMiilinl in lliit (lapri Of loiirtr, 
analyM* of any data ict tr<|inrrt, auovr all, 
ihr pioprr liaininK of uirnlihc qnr%lHini. 
Aim, ihr ri|(hi kind* of data iiuiu lir a<(r»»i- 
Mr lo aiuwrr ihr i|iirtiiont. llowr*rr, wr do 
frri ihal nilrrailivr di)(ilal iniagr pMMr*«ing 
ol a *aiiriy of groyiraphNally orirnirti data 
irit ran trad lu tub>iaiilial advanir* in oui 
iiiidrrnandniK of ihr catih 

Acknuwlrd|(mrntk 

I In* wuik waicaiiird oul with InnilniK 
Ituin llir ,\AS\ Olhir of S|UKr Sniriur anil 
A|iplHaiion*, NuniriirwaMr fa'ih Kr%oiiiir*, 
llnou|{li a conirail ('.I.S.S9.N9) Ironi ihr |cl 
Piii|nil*ion l,alMiraloiy. I’liHCtklliK w.it ion- 
diKird wiih a l*l)l‘-l |/.'I4 *yurin liMatrd wiili- 
III ihr \Va*lini)(ion I nivrrniy Kr^ional Planr- 
lary lniaf(r faiililv. I li.ink* air rxirndi'd lo 
Si. I.oiiit I'niiriiily lot *ri*niir il.iia, lo ihr 
Gropo*ilinnal Dnmoii, Drirnm .Ma|i|iiiiK 
ArrmiMir .A|{rn<.y (arnirr. Si. laniit. Mo., loi 
Ktaviiy d.iia, and In NO \.A for ihr io|H>)(ia- 
ph\ dala. 
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• rtcurd* on 
tompultr compalibir 
nvagnriic lapr, 
ANSI/ECMA formal 




• firld provrn 
in application* likr rnergy, 
trantporulion, agriculturr and 
environmental revearch 


• fully portable, operate* 
on rechargeable ballerie* 


• up lo 
too channel* 
f analog and/ 
or digital inputs 
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Oodata’* MI600L i* believed to be the only data logger on the market '1 
today with completely self-contained, data recording and replay capabili- j 
lie*. For scientiit* and engineer*, thi* mean*: J 


Convenient, fiilly-portable operation • MI600L operate* on it* own inter- 'j 

nal batteries with more than sufficient capacity to fully retvKd and replay *• 

each cartridge between recharging. The unit i* housed in a rugged, j 
weather resistant case. ] 


Accurate, versatile performance • MI600L uses individual conditioning 
cards, for each of 20 input channels (c.xpandable up to 100). The user 
configures lo his specific requirement. When needs change, cards are 
available at low cost to reconfiguie the system. 


Fast, convenient data processing • M 1600L records data in the ANSI/ECMA 
format on sMndard '/i-inch tape which can be replayed at high speed 
directly into a computer. Fur on-site verification the logger has built-in 
replay and display capabilities. 


For specifications, price and delivery, write or call; 
Crodata Corp. 

2SS Bear Hill Ru.id 
Waltham, MA 02154 
Tel: 617-890-0145 
Telex: 951693 1 
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